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Introduction and Motivation 
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• Structural (S uh . & HCF) assessment critical for turbomachinery 
flow path components undergoing possible resonance. 

• Resonance generally avoided, but impossible for higher modes 
found with modern analysis, especially with wide speed 
ranges. 

- Space Launch System upper stage J2-X Lox-H2 Engine Fuel Pump 
turbine stator operates from 26Krpm-34Krpm; 69N forcing excites 
modes 10-18 between 30KHz-40Khz. 

• Criteria triggers forced response analysis at worst case 
resonant condition. 

• Finite life analysis, where actual fatigue damage during 
operational time is calculated, frequently used if endurance 
limit criteria violated. 


Many Turbopumps "Dither" 
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May be beneficial to incorporate fact that real turbopumps 
dither about a nominal mean speed. 



J2-X Powerpack 
Adjusted Speed 
Trace 


During time speed is not exactly at natural frequency, damage 
accumulation is significantly reduced. 



Agenda 
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1. Introduction 

a) Motivation 

b) Literature Survey, Purpose Statement 

c) Characteristics of Excitation Frequency 

2. Theory and Development of Analysis Methods for Measured 
Speed Time Histories. 

a) Numerical Method, Time Step Convergence Study 

b) Analytical Method 

c) Calculation of Dither Life Ratio 

3. Monte Carlo Method for Unknown Speed History (Design), 

a) Sensitivity to Damping, Speed standard deviation 

4. Conclusion 


Literature, Purpose 
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Intial studies of response of systems with time varying excitation frequency Q 
by Lewis- 1932, Cronin- 1965. 

Lollack, 2002, defined reduction in peak response for monotically varying Q, 
useful for defining rate of sine-sweep tests. 

Henson, 2008, studied harmonically varying Q. 

For rocket engines, Q varies non-deterministically. Motivated previous work by 
authors (2010) that developed numerical approach for calculating response and 
general sensitivities. 

Unacceptable HCF factor for J2-X stator resonant 30Khz mode prompted need 
for practical technique. 


Purpose of this research 

- to develop procticol design techniques thot occountfor 
excitotion frequency stochasticity in the fatigue life of 
turbomachinery components. 


Excitation Frequency (kHz) 



Excitation Data 
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• Taken from hot-fire testing of J2-X and SSME. 

• Q = engine speed (hz)* [fo rcing pressure distortions/Rev] (FPR). 

• Since purpose is to examine fatigue life at resonance, actual mean speed 
adjusted to natural frequency for analysis. 

• Histograms for two different engines show ~ Gaussian distribution of 
speed. 




Coefficent of Variation (%) 


Theoretical Basis, Numerical Transient Solution 
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SDOF EoM 


x + 2 gcox + ry 2 x = 


m 


where 


m 


Q is derivative of 4>(t), constant in classical vibration analysis. 

For specified time-varying fl, t 

m = | Q(r)Jr 

o 

Calculate A necessary to generate peak resonant value of a a , t previously 
obtained by FEA, 

A 


°alt X 2 


co z 2£ 


• Now can solve for a a!t in EOM with using numerical Runge-Kutte procedure 
implemented in Matlab; agrees with Lollack's results for linearly varying Q. 




Finally, Calculate damage fraction O using Miner's rule, O = 
becomes 




dr 


, which 


Peak Dynamic Stress (psi) 



Convergence of Time Step in Transient Solution 
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• Previous work indicated convergence at At=l/40f n . 

• Initial studies here showed high freqency oscillation, so response 
and damage convergence studies performed -> At=l/120f n . 


104 Numerical Transient Response 



Time (s) 


• kt 2 Damage Accumulation during Test 



Time (s) 


Analytical Solution 
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Hypothesis from previous work that if /„□ 


d(speed) 
dt 


, then 


closed-form (computationally fast) standard analytical equation 
for SDOF steady-state response would be accurate. 


X 


CO 


steady-state 




Validation by comparing 
response with numerical 
solution. 


• 1CT 


Stress Response Zoom-In 



Analytical Solution Validation 
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Validation also shown in damage accumulation plot; error in 
analytical steady-state method is <1% (At < l/120f n required). 



Time (s) 


FFT of Speed also shows Analytical Sol'n Validity 
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. .... r dispeed) 

This assumption good with high FPR, driving /»/ — j t — ratio up. 

FFT of speed shows mostly below 100 hz, very low compared 
with natural frequency. 



"Dither Life Ratio" for Specified Excitation History 
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Calculation of damage performed considering dither for specific 
10 sec. window. 

Damage calculation assuming constant resonant excitation -> 
2.135 times more damage, call it "Dither Life Ratio". 


•1(T 2 




Monte Carlo for Unknown Frequency History 
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• During design phase, actual speed time histories unknown, but 
statistics from similar engines known. 

• Prompted development of Monte Carlo method using rapid 
analytical solution. 

• Speed vector created using 
Normal statistical distribution. 

• Powerpack data -> std dev 

=38.6 hz (cov=0.129%). = 

• MC results linear because rate | 

of change of frequency § 

variation not correct (and very "o 

high), but damage 
accumulation is accurate on 
the average. 

Time (s) 





Sensitivity of DLR to speed COV and C, 
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• Larger for high COV for speed, since more time spent off-resonance. 

• Larger for small < since peaks are sharper and time spent off- 
resonance will have less response. 



Conclusions 
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• Numerical and Analytical methods developed to determine damage 
accumulation in specific engine components when speed variation 
included. 

• Dither Life Ratio shown to be well over factor of 2 for specific case. 

• Steady-State assumption shown to be accurate for most turbopump 
cases, allowing rapid calculation of DLR. 

• If hot-fire speed data unknown, Monte Carlo method developed that 
uses speed statistics for similar engines. 

• Application of techniques allow analyst to reduce both uncertainty 
and excess conservatism. 

• High values of DLR could allow previously unacceptable part to pass 
HCF criteria without redesign. 

• Given benefit and ease of implementation, recommend that any finite 
life turbomachine component analysis adopt these techniques. 


